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Abstract We consider reduction of dimension for non¬ 
linear dynamical systems. We demonstrate that in some 
cases, one can reduce a nonlinear system of equations 
into a single equation for one of the state variables, and 
this can be useful for computing the solution when us¬ 
ing a variety of analytical approaches. In the case where 
this reduction is possible, we employ differential elimi¬ 
nation to obtain the reduced system. While analytical, 
the approach is algorithmic, and is implemented in sym¬ 
bolic software such as MAPLE or SageMath. In other 
cases, the reduction cannot be performed strictly in 
terms of differential operators, and one obtains integro- 
differential operators, which may still be useful. In ei¬ 
ther case, one can use the reduced equation to both 
approximate solutions for the state variables and per¬ 
form chaos diagnostics more efficiently than could be 
done for the original higher-dimensional system, as well 
as to construct Lyapunov functions which help in the 
large-time study of the state variables. A number of 
chaotic and hyperchaotic dynamical systems are used 
as examples in order to motivate the approach. 

Keywords reduction of dimension • differential 
elimination • nonlinear dynamics • chaotic attractors • 
computation of chaos 
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1 Introduction 

Nonlinear dynamical systems are ubiquitous in mathe¬ 
matics, engineering, and the sciences, with many real- 
world phenomenon governed by such nonlinear processes. 
In particular, non-equilibrium and chaotic dynamics are 
a continuing area of active research for applied mathe¬ 
maticians, as approximating such dynamics accurately 
and efficiently can be quite challenging. In the present 
paper, we shall consider reduction of dimension for non¬ 
linear dynamical systems. This approach has previously 
been employed in the literature in order to enable the 
construction of Lyapunov functions [1] and equilibrium 
dynamics [5], as well as to allow one to more easily 
approximate chaotic attractors analytically [31lll[51[n]. 
One method for reduction of dimension is differential 
elimination, in which one algorithmically reduces the 
nonlinear dynamical system into a single ordinary dif¬ 
ferential equation (ODE) for one of the state variables. 
However, this is possible only when the system reduces 
to an ODE; if the reduction is instead to an integro- 
differential equation, the process is not algorithmic and 
specific cases must be handled with more individual 
care. Our focus shall be on dynamical systems giv¬ 
ing chaotic dynamics, but the approach can certainly 
be applied for non-chaotic ODE systems. We give an 
overview of reduction of dimension, after which we demon¬ 
strate in several ways why one might wish to apply this 
technique. 

With the wide range of numerical methods available 
for solving nonlinear first-order ODE systems of even 
high order, one may wonder why it might be advanta¬ 
geous to convert such systems into a single higher-order 
ODE. We shall mention several situations in which the 
differential elimination, and more generally reduction 
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of dimension, may prove useful. We then outline the 
paper. 

Often times, if one is trying to approximate the so¬ 
lution to a nonlinear system through some sort of ana¬ 
lytical approximation, via series, perturbation, or more 
complicated approaches, one quickly finds that the cou¬ 
pled equations require balancing many terms coming 
from the expansion for each of the state variables. In the 
case of a single state variable governed by a higher-order 
ODE, one need only track terms in a single asymptotic 
expansion. This approach has been applied when using 
Taylor series, approximate Fourier series, and asymp¬ 
totic expansions in other types of basis functions to 
the solution of a system of nonlinear ODE. In hybrid 
analytic-numeric methods, such as the homotopy anal¬ 
ysis method such reductions of a system to a 

single equation also simplifies the optimization problem 
which is solved to obtain the error-minimizing solution 
(see, for instance, [5], where the present approach is 
used in such a capacity). Therefore, the reduction of 
order can greatly reduce the complexity of analytical 
calculations under several frameworks. 

Contraction maps or Lyapunov functions are useful 
tools for discussing the convergence of solutions to non¬ 
linear dynamical systems to large-time steady or quasi¬ 
steady dynamics. In situations where contraction maps 
or Lyapunov functions are known for a given dynamical 
system, the state variable governed by a single higher- 
order ODE necessarily results in a contraction map in 
the single state variable. However, as is well known to 
those studying stability of nonlinear systems, it is not 
often easy to obtain contraction maps for complicated 
systems. As we shall show here, it is possible to use 
the reduction of a system to a single higher-order ODE 
in order to construct a contraction map for the state 
variable governed by the aforementioned higher-order 
ODE. The existence of such a map can then be used to 
deduce the large-time dynamics of the state variable, 
as well as for the other state variables in the original 
system. One example of this is given in [1], and other 
examples are provided in Section 5. 

Related to both the topic of analytical approxima¬ 
tions and Lyapunov functions would be long time dy¬ 
namics and equilibrium behavior of nonlinear dynami¬ 
cal systems. Indeed, in order to study the equilibrium 
structure of a high-order system of ODEs, one must 
solve a coupled system of nonlinear algebraic equations 
in order to recover the fixed points for the state vari¬ 
ables. First reducing the system to a single ODE allows 
one to obtain a single nonlinear algebraic equation for 
the fixed point of a single state variable, which can then 
be used to recover the fixed points of the other state 
variables. Therefore, when such a reduction to a single 


ODE is possible, the need to solve a nonlinear alge¬ 
braic system for all of the fixed points simultaneously 
is eliminated, resulting in what is often a far less com¬ 
putationally demanding problem. 

Another topic is great recent interest in nonlinear 
science has been both the synchronization of chaos nni 
[nni] and the control of chaos. In situations where one 
is interested in mitigating the possibility of emergent 
chaos, one can couple a chaotic system to various con¬ 
trol terms, or indeed to additional dynamical systems, 
which may lend a degree to stability. Under such ap¬ 
proaches, one often increases the complexity or even 
the dimension of the dynamical system being solved. 
As such, methods to reduce the dimension of such sys¬ 
tems could improve compatibility. Furthermore, since 
the control of chaos is often linked to a control term 
which itself is determined by a Lyapunov function, the 
construction of contraction maps through the reduction 
approach outlined here could be of great use. 

As stated before na, the competitive modes anal¬ 
ysis gives an interesting link between the geometry of 
phase space possibly yielding chaotic trajectories (recall 
that the competitive modes requirements appear to be 
a necessary, albeit not sufficient, condition for chaos 
[m[TSl[TBlll7lll8l[Ta i. Conversely, the differential elimi¬ 
nation may cast light into the geometry of solutions in 
the space of derivations. Since this result of the differen¬ 
tial elimination is a single higher-order ODE, and since 
any chaos emergent from the nonlinear system should 
be encoded in the single higher-order ODE, the differ¬ 
ential algebraic structure of such equations may cast 
light into practical geometric tools by which one may 
study systems in which chaos is observed. In particu¬ 
lar, through this reduction approach, the calculation of 
mode frequencies in the standard competitive modes 
analysis becomes much simpler. 

This paper is outlined as follows. In Section 2, we 
provide an algorithmic approach, to differential elim¬ 
ination for nonlinear dynamical systems based upon 
differential algebra. First laying out the general the¬ 
ory, we then give specific MAPLE code for performing 
the differential elimination in a systematic manner. The 
algorithmic approach is useful in the case where the 
dynamical system can be reduced to a single ODE in 
terms of only one of the state variables. In Section 3 we 
implement the approach in order to reduce a variety 
of chaotic and hyperchaotic systems, finding that the 
form of the nonlinearity in the dynamical system will 
strongly influence the reducibility properties. However, 
in cases where the dynamical system is not reducible us¬ 
ing differential elimination, one may still obtain more 
complicated reductions, for instance in terms of inte¬ 
grals, resulting in more complicated integro-differential 
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equations for the reduced state variable. The possible 
results are illustrated through concrete examples for 
the Rossler system (which is completely reducible), the 
Lorenz system (which is partially reducible - that is, 
reducible in some but not all state variables), and the 
Qi-Chen-Du-Chen-Yuan (which is irreducible under dif¬ 
ferential elimination, but which can be reduced to an 
integro-differential equation). We give summarizing ob¬ 
servations regarding the reducibility of dynamical sys¬ 
tems in Section 4. 

The remainder of the paper is devoted to applica¬ 
tions of reduction of dimension for dynamical systems. 
In Section 5, we demonstrate that reduction of dimen¬ 
sion can be useful for obtaining contraction maps and 
Lyapunov functions, which in turn may be used to de¬ 
termine asymptotic stability of dynamical systems and 
also to control chaos in such systems. In Section 6 we 
demonstrate that reduction of dimension can be used 
to simplify calculations involved in certain techniques 
for studying the solutions of nonlinear dynamical sys¬ 
tems. Indeed, when applicable, we find that the ap¬ 
proach greatly reduces the number of nonlinear alge¬ 
braic equations required to be solved when constructing 
trajectories in state space via undetermined coefficient 
methods by a factor of 1/n, where n is the dimension 
of the dynamical system, meaning that the number of 
equations needing to be solved will not increase with 
the size of the system. Furthermore, when applying the 
competitive modes analysis (which is a type of diagnos¬ 
tic criteria for finding chaotic trajectories in nonlinear 
dynamical systems), we find that only one binary com¬ 
parison is needed of one first reduces the dimension of 
the dynamical system so that there is a single equation 
for one state variable. In contrast, there are normally 
of order 2”“^ comparisons needed for an n-dimensional 
dynamical system. In Section 7 we provide a discussion 
and possible avenues for future work. 


2 Algebraic approach to differential elimination 

Systems of differential equations are ubiquitous and 
widely studied. Ritt [50] and Kolchin m pioneered 
the field of differential algebra, an algebraic theory for 
studying solutions of ordinary and partial differential 
equations. We are particularly interested in differential 
elimination, an algorithmic subtheory that can simplify 
systems of parameterized algebraic differential equa¬ 
tions. This permits one to reduce the dimension of a 
dynamical system so that one is left with a single ODE 
in the state variable. 


2.1 Algebra preliminaries 

Here we briefly review concepts from algebra and dif¬ 
ferential algebra. For reference books in differential al¬ 
gebra, see nnum. If / is a subset of a ring R then (!) is 
the (algebraic) ideal generated by I. Let I be an ideal of 
R. Then \/7 denotes the radical of I. A derivation over 
a ring i? is a map R^ R which satisfies (we write a is 
the derivative of a), for every a,b & R, (a + b) = d + b 
and (a b) = {d)b + a{b). The held of differential alge¬ 
bra is based on the concept of a differential ring (resp. 
held), which is a ring (resp. held) R endowed with a set 
of derivations that commutes pairwise. A differential 
ideal [/] of a differential ring R is an ideal of R stable 
under the action of derivation. 

Diherential algebra is more similar to commutative 
algebra than analysis. In commutative algebra, Buch- 
berger solved the membership problem (tests whether a 
given polynomial is contained in a given ideal) through 
the theory of Grobner bases [3^ . From algebraic geom¬ 
etry, we know the set of polynomials which vanish over 
the solutions of a given polynomial system form an ideal 
and even a radical ideal [31]. In the case of differential 
equations, the set of differential polynomials which van¬ 
ish over the analytic solutions of a given system of dif¬ 
ferential polynomial equations form a differential ideal 
and even a radical differential ideal [50] . Ritt solved the 
theoretical problem (of membership for radical differ¬ 
ential ideals) and developed algorithmic tools to solve 
systems of polynomial ODE and PDEs; however, Ritt’s 
algorithm requires factorization. 

Due to the complexity of factorization, Boulier and 
co-authors avoided it by developing the Rosenfeld-Grobner 
algorithm, based on the work of Seidenberg and Rosen¬ 
feld, and incorporating Grobner bases mMm- Since 
then, the algorithm has been improved both theoreti¬ 
cally and practically P1I551I5H] and it no longer requires 
Grobner bases. It is available in the Dif f erentialAlgebra 
package in MAPLE [55] and SageMath as an interface 
for the BLAD and BMI libraries mm- 

Algorithmically, differential elimination involves ma¬ 
nipulation of finite subsets of a differential polynomial 
ring R — K{U} where K is the differential field of coef- 
hcients (i.e. K = Q), and 17 is a finite set of dependent 
variables. The elements of R are differential polynomi¬ 
als, which are polynomials built over the infinite set 
of all derivatives 0U, of the dependent variables. Gon- 
sider a system S of polynomial differential equations, 
here, we consider the Lorenz system of three ordinary 
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differential equations: 
x'l = a{x2 — xi), 

X2 = Xi{b - X3) - X2 , ( 1 ) 

X 3 = X 1 X 2 — CX 3 . 

The Lorenz system can be re-written as: 

{ -x'l -I- a(x 2 - X 3 ) = 0 , 

-X 2 -b Xi(& - X 3 ) - X 2 = 0 , (2) 

—X 3 -b X 1 X 2 — CX 3 = 0 . 

The Rosenfeld-Grobner algorithm takes as an input a 
differential system E and a ranking. A ranking > is any 
total ordering over the set OU of all derivatives of the el¬ 
ements of U which satisfies the following axioms: a < a 
and a < & =b d < & for all a, 6 G OU. The Rosenfeld- 
Grobner algorithm transforms E into finitely many sys¬ 
tems called regular differential systems, which reduces 
differential problems to purely algebraic ones that are 
triangular. The next step is purely algebraic and trans¬ 
forms the regular differential system into finitely many 
characteristic presentations, Ci,.. .Cr. Rosenfeld-Grobner 
outputs this finite family Ci,.. .Cr of finite subsets of 
K{U}\K, where each Ci defines a differential ideal [Ci]. 
The radical \/[E] of the differential ideal generated by 
E is the intersection presented by characteristic sets: 

= [Q] n • • • n [Cr\. 

Note differential ideals [Ci] do not need to be prime, 
however by hazard’s lemma, they are necessarily rad¬ 
ical. Differential algebra elimination has proven useful 
for parameter estimation, identifiability, and model re¬ 
duction of biological and chemical systems [291128] . 

2.2 Gomputational Method 

We demonstrate reduction in dimension via differen¬ 
tial elimination algorithm Rosenf eldGroebner in the 
DifferentialAlgebrapackage implemented in MAPLE. 
For sake of using a concrete example, we choose the 
Lorenz system. First, we call the package: 
with(DifferentialAlgebra): 

Next we input the Lorenz system: 

sys := [-(diff(xl(t),t))+a*(x 2 (t)-xl(t)), 

-(diff(x2(t),t))+xl(t)*(b-x3(t))-x2(t) , 

-(diff(x3(t),t))+xl(t)*x2(t)-c*x3(t)] 

Next we form our differential ring, embedding the rank 
of dependent variables in blocks and independent vari¬ 
ables in derivations. Since we are considering ordi¬ 
nary differential equations, derivations is set to one or¬ 
dering, time t. We remark that the DifferentialAlgebra 


package enables differential elimination of PDEs by in¬ 
cluding additional inputs for the derivations (e.g., 
derivations= [u,x,t] . Note, sys is assumed to have 
coefficients in the field Q[xi, X2, X3] obtained by adjoin¬ 
ing the independent variables to the field of rationals 
and symbolic parameters a, b, c are considered arbitrary 
in the coefficient field. To form the differential ring, we 
input: 

R := DifferentialRingCblocks = [x3, x2, xl] , 
derivations = [t]) 

Note that xi stands to the rightmost place on the list 
which identifies that we are attempting to reduce the 
differential equation to only one variable, i.e. xi(t). This 
ranking eliminates X 3 with respect to X 2 , and then elim¬ 
inates X 2 with respect to xi. We now call the Rosenfeld 
Grobner algorithm for our system and differential poly¬ 
nomial ring: 

G := RosenfeldGroebner(sys, R) 
simplify(Equations(G[l], solved)) 

This will return the characteristic presentation (which 
should be understood as an intersection), with the equa¬ 
tions given by the ranking, with the final equation a 
single ODE for xl(t), provided that it exists and can 
be computed by the algorithm. In some cases the algo¬ 
rithm will keep running and therefore should eventually 
be terminated by the user. For such cases, it is unlikely 
that a reduction of the specified form exists. However, 
as we shall consider in the next section, when the reduc¬ 
tion is to an integro-differential equation, rather than 
an ODE, the approach will not identify the reduced 
equation. 

3 Reduction of Dimension: Applications 

Here we apply the method of differential elimination 
to several nonlinear dynamical systems known to give 
chaos, in order to see if these equations can be reduced. 
We hrst apply the algorithmic approach outlined in 
Section 2, finding that the approach gives a complete 
reduction (all state variables can be isolated and ex¬ 
pressed as the solution to single uncoupled ODEs), a 
partial reduction (one or more, but not all, state vari¬ 
ables can be isolated and expressed as the solution to 
single uncoupled ODEs), or returns no reduction (the 
algorithm does not complete in a fixed amount of time), 
in which case none of the state variables can be ex¬ 
pressed as a solution to a single ODE reducible from 
the original system. For simplicity, we shall only con¬ 
sider autonomous systems. 

We consider a number of examples of chaotic sys¬ 
tems in Table 1, with the results of the differential elim¬ 
ination algorithm given. We also give a summary of the 
dynamics of the example equations selected. Since the 
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form of these equations may vary through the literature, 
we give a list of the specific form of the equations con¬ 
sidered, in Appendix A. Note that we have considered 
the differential elimination algorithm for the arbitrary 
parameter values listed in Appendix A. Table 1 demon¬ 
strates that the structure of the dynamical system tends 
to play a strong role in whether the system can be re¬ 
duced. Indeed, equations with a single nonlinearity tend 
to be completely or partially reducible, hence at least 
one state variable can be solved for via a single nonlin¬ 
ear ODE. On the other hand, the equations with many 
nonlinear terms or higher-order degree of nonlinearity 
(we consider only equations with polynomial nonlinear¬ 
ities) tend more often to be irreducible using the ap¬ 
proach. Of the listed equations, note that the Rossler 
system is one of the few completely reducible systems, 
lending validity to the belief that it is indeed one of the 
simplest possible continuous- time dynamical systems 
giving chaos. Meanwhile, the commonly studied Lorenz 
system is only partially reducible under the approach. 
More complicated systems tend to be irreducible under 
the algorithm, and many of these give more complicated 
dynamics such as multiple scroll attractors. 

Note that the algorithm returns a ‘No’ if a reduction 
is not obtained within a given time interval. For cases 
where the algorithm found a reduction, the computa¬ 
tion time was fairly quick. We are therefore comfortable 
in assuming that a reduction to an ODE does not exist 
in cases where the the algorithm times out. For such 
cases, the system may still admit a reduction, but not 
strictly in derivatives of one of the state variables. One 
such example would be a system which is reducible to 
an integro-differential equation in one of the state vari¬ 
ables, but never to simply an ODE. 

We next consider hyperchaotic systems (chaotic sys¬ 
tems giving two or more positive Lyapunov exponents) 
in Table 2. Again, we find that the more complicated 
the functional form of the nonlinarities, the less likely 
a system seems to be reducible. Furthermore, hyper¬ 
chaotic generalizations of known chaotic systems ap¬ 
pear to maintain their reducibility properties, since of¬ 
ten a simple additional equation is added to make a 
chaotic system hyperchaotic. The hyperchaotic Rossler 
system is completely reducible, as was the related chaotic 
Rossler system, again suggesting that the chaotic and 
hyperchaotic Rossler systems are some of the simplest 
systems which still exhibit chaos and hyperchaos, re¬ 
spectively. 

The results indicate that completely reducible sys¬ 
tems are perhaps the simplest systems giving chaos 
or hyperchaos. Again, this would support the qualita¬ 
tive and topological claims that the Rossler systems 
are some of the simplest possible equations permitting 


System 

Dynamics 

Reducible? 

Lorenzl34 11351 

3D: 1-2-2 

Partial 

Modified Chua’s circuit I36II371I3SI 

3D:3-1-1 

Gomplete 

Chen-1 ;eel39ll4f)l 

3D:2-2-2 

Partial 

Rabinovich-Fabrikant I4l 11421143 1 

3D:3-3-3 

No 

R6sslerl44 11451 

3D:l-l-2 

Gomplete 

Chenl46 11471 

3D: 1-2-2 

Partial 

T, 51480491 

3D: 1-2-2 

Partial 

T-system 15905111521 

3D: 1-2-2 

Partial 

Qi-Du-Chen-Chen-YuanlF),'II 

4D:3-3-3-3 

No 

Qi-Chen-Du-Chen-Yuanl54l 

3D:2-2-2 

No 

Generalized Lorenz|55||56||57| 

3D:2-2-2 

No 

Blue-sky catastroDhel580590fi0l 

3D: 3-3-3 

No 

Lorenz-Stenflolfil 11620631 

4D: 1-2-2-1 

Partial 

Genesio-Tesi I64U65I 

3D:l-l-2 

Partial 

Arneodo-Coullet-Tresser[66J 

3D:l-l-3 

Partial 


Table 1 List of chaotic systems and their reduction prop¬ 
erties. The numbers in the dynamics column indicate the 
dimension of the system and then degree of each poly¬ 
nomial in the respective reaction functions. For instance, 
if i = A{x,y,z), y = B{x,y,z), z = C{x,y,z), then 
3D:deg(A) — deg(i?) — deglC) is reported, where deg(A) de¬ 
notes the degree of A, and so on. When the system has the 
property that it may be reduced to a single ODE in any state 
variable, we say that it is completely reducible, and record 
a ‘Complete’. If the system may be reduced to an ODE in 
one or more, but not all, state variables, we say the systems 
is partially reducible, and record a ‘Partial’. Finally, when a 
system is not reducible to a single ODE in any state vari¬ 
able, we record a ‘No’. We note that specific forms of some 
equations can change from paper to paper, so we record the 
specific equations used in Appendix A. 


System 

Dynamics 

Reducible? 

Rossler 1671 

4D:l-l-2-l 

Complete 

Chen 1680691 

4D:l-2-2-l 

Partial 

Lii CO] 

4D: 1-2-2-2 

No 

Modified Lii 1711 

4D:2-2-2-l 

No 

Wang-Liu [72J 

4D:l-2-2-l 

Partial 

,lia 1730741 

4D: 1-2-2-2 

Partial 

QWWC system PtMTHI 

4D:2-2-2-2 

No 


Table 2 List of hyperchaotic systems and their reduction 
properties. The labeling is the same as was given in Table 
1. We note that specific forms of some equations can change 
from paper to paper, so we record the specific equations used 
in Appendix B. 


chaos m, as they each involve only a single quadratic 
nonlinearity. On the other hand, systems with stronger 
polynomial nonlinearities, or systems with many non¬ 
linear terms, appear to often be irreducible under dif¬ 
ferential elimination. Note that for cases where the re¬ 
duction might involve integrals, resulting in a type of 
integro-differential equation, the differential elimination 
algorithm would miss such a reduction, even though it 
exists. This is due to the fact that the differential elim¬ 
ination algorithm is working over the ring of deriva¬ 
tions, which does not include integrals. Indeed, since 
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integral operators are fairly cumbersome to introduce 
compared to their differential operator counterparts (we 
discuss this later in Section 7), obtaining an algorith¬ 
mic approach including integrals would be challenging. 
Therefore, the differential elimination algorithm out¬ 
lined in Section 2 appears to be a very useful tool for 
reducing the dimension of dynamical systems, provided 
that a reduction to a single ODE exists. For the more 
complicated models, we find the need to proceed on a 
case-by-case basis with manual manipulations due to 
any integration needed. 

We demonstrate reduction of dimension for chaotic 
systems into single higher-order ODEs in the next three 
subsections. We pick a case where all state variables can 
be isolated (the Rossler system), a case where one of the 
state variables can be isolated in terms of a differential 
equation (the Lorenz system), and finally a case where 
none of the state variables can be isolated in terms of a 
differential equation (the Qi-Chen-Du-Chen-Yuan sys¬ 
tem) so that any reduction would necessarily involve in¬ 
tegrals. For all cases considered, we let x,y,z £ C'"'(K) 
where n is the dimension of the relevant dynamical sys¬ 
tem, and we take a, 6, c G M to be parameters. 


3.1 Rossler System 

The Rossler equation [^H5] reads 

X = -y- z, 

y = x + ay, (3) 

z = b + x(x — c). 

We first obtain the ODE for y{t). Note from the 
second equation that x = y —ay, so that y = y —ay and 
hence from the first equation we have z = —y + ay — y. 
Placing these into the third equation, and performing 
algebraic manipulations, we obtain 

y - ay+ y - {y - ay + y) {y - ay - c)+ b = 0 . (4) 

Note that this equation is third order, and therefore the 
information of the three-dimensional system m can be 
encoded in this single ODE. By similar manipulations, 
one may arrive at an equation for x{t), 


{a + c — x)^ 


i- 

\dt 


(x 



X — ax + X + b 


a + c — X 



= 0 , 
( 5 ) 


and an equation for z{t), 

z^ (— a^ + c") -—- + z^z — az^ + cz^ = 0 . (6) 

at / z 


3.2 Lorenz System 


The Lorenz system mM is given by 
x = a{y-x), 

y = x{b-z)-y, (7) 

z = xy — cz . 


Observing from the first two equations that 
1 

y = X -\—X 
a 

and 


( 8 ) 


, x + {l + a)x + x 

z = h -, (9) 

ax 

the third equation can be used to obtain a single ODE 
for the state variable x{t), viz.. 



X -b (1 -|- a)x -b X 

X 


-b ax'^ + x^x — abcx^ = 0 . 


( 10 ) 


This agrees with what one obtains from the differential 
elimination. On the other hand, we observe that the 
algorithmic approach to differential elimination is use¬ 
ful for situations in which there is no obvious route to 
reduce a system into a single equation (through elim¬ 
inations and substitutions). A good example of this is 
found when trying to obtain a differential equation for 
the state variable z{t) alone. Using the differential elim¬ 
ination, we arrive at a rather complicated equation of 
the form 


{b - z){'z)'^ + Pl{z, z, z)'z + P 2 {z, z, z) =0, (11) 

where Pi and P 2 are complicated polynomials that we 
do not list for sake of brevity. Interestingly, this is a fully 
nonlinear equation, since the highest order derivative 
enters into the equation nonlinearly. In contrast, the 
equation obtained for the state variable x{t) is quasi- 
linear, since it is linear in the highest derivative. One 
could differentiate the equation for z{t) in order to iso¬ 
late the highest derivative, but by doing so one would 
increase the differential order of the system, thereby 
decreasing the regularity of the system. This is partic¬ 
ularly important in cases where the solution z(t) may 
only be C^(]R). 

When a system is nonlinear, there may of course be 
forms of the nonlinearity which do not permit one to 
obtain an equation for a single state variable in terms of 
that state variable and its derivatives. A good example 
of this is the state variable y{t) in the Lorenz system. 
The algorithmic differential elimination finds no closed 
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differential equation for y(t). As it turns out, the rea¬ 
son for this is that any equation governing y{t) alone 
will necessarily involve integral terms which cannot be 
eliminated (due to the nonlinearity of the equation). To 
see this, note that if we consider the first equation in 
the Lorenz system, which may be written in the form 
{e°‘^x)' = ae°'*y, we find 


:(f) = xq + ae / e°'^y{s)ds . 

Jo 


( 12 ) 


Here Xq is the initial value of the state x{t), that is 
a;(0) = Xq- Yet, from the second equation in the Lorenz 
system, we have z = b — (y + y)lx, which yields 


z{t) = b — 


y + y 


Xq + /g e°‘^y{s)ds 


(13) 


Placing the representations for x{t) and z{t) into the 
third equation of the Lorenz system, and performing al¬ 
gebraic manipulations to simplify the resulting expres¬ 
sion, we obtain 

(y + (1 + c)y + cy) ^a;o -b ae““* e°^''y(s)d. 

+ {y + y) {ay - J e°‘’'y{s)d. 


y [ xo + ae 


— cb i Xq + ae 


e°'^y{s)ds 


(14) 


-j\-^y{s)ds^ = 0 . 


Note that the equation both involves an integral and is 
non-autonomous. 


3.3 Qi-Chen-Du-Chen-Yuan System 

We now consider the Qi-Chen-Du-Chen-Yuan (QCDCY) 
system [54) . which is given by 

X = a{y - x) + yz , 

y = bx — y — xz , (15) 

z = xy — cz . 

The differential elimination algorithm indicates there is 
no reduction to a single ODE in any of the three state 
variables. This system has a quadratic nonlinearity in 
each equation, and this added complication is behind 
the difficulties in obtaining such a reduction. However, 
we may still obtain an equation for a single state vari¬ 
able, if we are willing to include integral terms. Due to 
the complexity in obtaining such an equation, we shall 
restrict our attention to finding a single equation for 
the state variable z{t), noting that similar approaches 


can be used to hnd a single equation for either of the 
other two state variables, x{t) or y{t). 

Let us begin by noting that the second equation in 
the QCDCY system implies (e*?/)' = e*{b — z)x, while 
placing this into the third equation in the QCDCY sys¬ 
tem gives (a + z)(z + cz) = a;e““‘(e“*a;)' = x{x -b ax). 
This, in turn, implies that state variables x{t) and z(t) 
satisfy 


e^“*(a;(f))^ = Xq- b 2 f e‘^°''^{a + z{s)){z{s) + cz{s))ds . 

Jo 

(16) 

where x(0) = xq. The first equation in the QCDCY 
system has not been used, and we place this relation 
into that equation to obtain a single equation for the 
state variable z{t). After several algebraic and differ¬ 
ential manipulations, we arrive at the single equation 


2 e ^“‘ — a - (^0 ^0 

-b 2 (xq -b J[z, z]) ^ (e^“‘(i -b cz){a + z)) 

— 2e'^“‘(i -b cz)‘^{a + z)"^ 

— 2(6 — z)(a -b z) (xq -b J[z, i]s)^ = 0 , 

(17) 

where we have defined the integral operator 
J[z,i;]=2 f e^°'^{a +z{s)){z{s) + cz{s))ds . (18) 

Jo 

Similar results can be obtained for the other state 
variables. The fact that the obtained equations involve 
an integral operator which cannot simply be differen¬ 
tiated away demonstrates why the differential elimina¬ 
tion algorithm was not useful for this case. Still, per¬ 
forming the manipulations by hand, we have reduced 
the fairly complicated QCDCY system into a single 
integro-differential equation, thereby reducing the di¬ 
mension of the original system. 


4 Reductions of n-dimensional dynamical 
systems 

We now give some summarizing remarks based on what 
we have seen in the previous sections. We shall assume 
that each system is coupled through at least one state 
variable (otherwise the state variables naturally sepa¬ 
rate into distinct lower-dimension equations, and the 
approach is not needed). 
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4.1 Linear Systems 

For first order linear systems of dimension n, there is 
always a reduction into a single higher-order ODE. This 
follows from the process of Gaussian elimination. In the 
case that the matrix of coefficients for such a first-order 
system is full rank, the resulting higher-order ODE will 
be of order n. If the matrix of coefficients is singular, 
then the resulting higher-order ODE will be of order 
less than n. 


4.2 Reducible Nonlinear Systems 

For first order nonlinear systems of dimension n, there 
are multiple possibilities, owing to the structure of the 
nonlinearity. 

In cases where the system permits the complete 
differential elimination (an example being the Rossler 
equation), all state variables in a first order nonlin¬ 
ear system can be expressed in terms of a higher-order 
ODE. Note, however, that it is possible for the order 
of the single ODE to be different from the dimension n 
of the first order system. As an example of this point, 
consider the system 

X = X — y — z , 

y = x‘^, (19) 

Z = X — x^ . 

Clearly, differentiation of the first equation gives x = 
X — y — z = x — x^ + x — x^. So, we obtain 

X — X — x + x^+x^ = 0, (20) 

which is a second order equation for the state variable 
x{t), even though the original system was first order. A 
similar example can be found in [I] , where a fourth or¬ 
der nonlinear dynamical system was reduced to a single 
second order nonlinear ODE. 

It is possible for a system to be reduced to a single 
equation, which is not an ODE. This was evident even 
for the Lorenz equation, where an equation for one of 
the state variables involves an integral term in addition 
to derivative terms. Note that the equation was not 
closed under any number of differentiations, due to the 
form of the intergral terms. As such, the single reduced 
equation for the state variable could never be expressed 
strictly as an ODE of any finite order. Note that this 
can occur for one of the state variables, while a dif¬ 
ferent state variable might satisfy a finite order ODE. 
For such cases, the nonlinearity in the system results in 
their being certain favored state variables with which 
to perform the reduction to a single ODE. 


4.3 Differentially Irreducible Nonlinear Systems 

We have observed that for more complicated nonlin¬ 
ear dynamical systems, there is no reduction to a single 
ODE in one state variable. While it may be the case 
that differential elimination does not pick up an ODE 
that does exist, it seems as though the failure of dif¬ 
ferential elimination is a sign that integrations will be 
needed in order to reduce the dimension of such sys¬ 
tems. Indeed, when integrations of this kind are called 
for, the manipulations are no longer confined to the 
specified differential ring, and the differential elimina¬ 
tion cannot be performed. While one can attempt these 
integrations manually, as opposed to algorithmically, 
obviously it would be desirable to have some kind of 
algorithmic approach. Perhaps one may adjoin inte¬ 
gral operators to the differential ring, in order to per¬ 
form reductions for more complicated nonlinear sys¬ 
tems. This would likely work in cases like the Lorenz 
system, for which there is partial reducibility under dif¬ 
ferential elimination. Eor instance, if one were to define 
a new variable Y(t) = f* e°'^y(s)ds, then one would 
obtain a non-autonomous ODE for Y (t) from equation 
m- Therefore, this fairly simple integral transforma¬ 
tion, in addition to differential operators, can reduce 
the dimension of the Lorenz system with respect to 
the state variable y{t). However, in cases like that of 
the QCDCY system, note that the form of the integral 
operator given in (HID is rather complicated, depend¬ 
ing nonlinearly on both the state variable z{t) and its 
derivative z{t). For such cases, there is no combination 
of elementary integral transforms that can be adjoined 
to the differential ring which would permit reduction 
of dimension to a single ODE. As such, it appears as 
though reduction of dimension for certain more com¬ 
plicated systems will result in reductions to integro- 
differential equations, rather than ODEs, for some fun¬ 
damental reason related to how complicated the orig¬ 
inal dynamical system is. Therefore, the study of pos¬ 
sible algorithmic methods for the reduction of dimen¬ 
sion for dynamical systems into single scalar integro- 
differential equations would be an interesting and po¬ 
tentially very useful area of future work. 

5 Contraction Maps and Lyapunov Functions 

Turning out attention now to practical applications for 
reduction of dimension, recall that contraction maps 
and Lyapunov functions are useful tools for studying 
the asymptotic stability of nonlinear dynamical sys¬ 
tems. In this section, we use the three examples worked 
explicitly in Section 3 in order to demonstrate the util¬ 
ity of reduction of dimension for finding Lyapunov func- 
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tions. Using these results, we can recover stability re¬ 
sults for these dynamical systems which were obtained 
through other approaches, and which agree with exist¬ 
ing results in the literature. 


5.1 Rossler System 

The Rossler system has two equilibrium values, ±t/*, 
for t/(t), and the constant y* must satisfy the quadratic 
equation 

a{y*f + cy* + b = Q. ( 21 ) 

In order to discover a Lyapunov function for the 
Rossler system, it is tempting to assume a bowl-shaped 
map of the form ax^ -\- jSy'^ + 72 :^, or minor variations 
on this theme involving higher power polynomials of 
even order, but the approach evidently proves fruitless. 
Therefore, we shall use one of the three equations ob¬ 
tained for the isolated state variables of the Rossler 
system. 

Consider equation ® for the Rossler system @. 
Let us write Y{t) = y{t) — y* in the neighborhood of 
either equilibrium value y*. This transformation will 
prove useful, as the Lyapunov function needs to vanish 
at the equilibrium value selected. (There is therefore the 
need to construct such a function in a neighborhood of 
each equilibrium point.) Under this transformation, 0 
is put into the form 

Y - aY + Y - (y - aY + Y^ (y - aY^ = 0 . (22) 

Let us define a function m = Y — oY -|- U so that (|22|1 
is put into the form 

m-{Y - aY)m = 0 . (23) 

Observe that (ESI can be written as 
m-e“‘(e"“‘U)'m = 0. (24) 

From this, we recover 

Y — aY + Y = m = moexp e°‘’(e~°‘'’Y(()yd('j , 

(25) 

where toq is a constant of integration. As we are inter¬ 
ested in recovering information about the asymptotic 
stability of the Rossler system, let us pick the initial 
condition Y (0) = e. This corresponds to setting the ini¬ 
tial condition such that it is contained within a neigh¬ 
borhood of the equilibrium value. Let us also restrict 


|a| < 2 (this will simplify the mathematics, and is con¬ 
sistent with the physics of the Rossler system). Then, 
we obtain 


Y(t) = 


— ppat/2 


ee ' < cos 


\/4 — 


t\ -f C sin 


'JA — a? 


-I-Too J K{t,s)exp e“^(e “^F(C))'dC^ ds , 

(26) 


where C is a constant that will depend on the initial 
value of Y ( 0 ) (the value of which will not impact our 
analysis) and K{t,s) is the kernel 

A:(t,s) = sin s)^ . (27) 

Observe that for — 2 < a < 0 the map is a contraction. 
Given arbitrarily small e > 0, for large enough time 
t(e) > 0, the solution Y{t) will lie in a neighborhood 
—e<Y{t)<e for all t > f(e). Therefore, F —>■ 0 as 
t —>■ 00. Yet, by definition of Y{t), this implies y —>■ 
y* as t —>■ 00. Using this, one may shown that x —>■ 
—ay* and z —>■ —y* as t — ^ 00. Hence, we have shown 
that a < 0 gives a stable solution, which was already 
known from different work. The nice thing about this 
approach is that is allows us to bypass a linear stability 
analysis involving the calculation of eigenvalues at the 
algebraic solution to y* found from m- Indeed, we 
did not even need to calculate the equilibrium value y* 
for the present analysis, as the analysis holds for an 
arbitrary equilibrium value satisfying (EH). 


5.2 Lorenz System 

In order to find a Lyapunov function for the Lorenz sys¬ 
tem in a neighborhood of the zero equilibrium {x, y, z) = 
(0, 0, 0), let us assume a bowl type function of the form 

V{x,y,z) = ax'^ + Py^^+-fz'^ , (28) 

where a > 0, /3 > 0, and 7 > 0 are constant param¬ 
eters to be selected, recall that physically interesting 
model parameters a, b, and c are positive. Then, the 
time derivative of V is given by 

-V = —aax'^ — Py"^ — 'ycz'^ + {aa + pb)xy + {'y — Pb)xyz . 

(29) 


Clearly, we should take 7 = pb. Note that 
— {y/aax — \/PyY = —aax^ — Py'^ + 2y/apaxy . (30) 
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Then, 

-V = —{\/7mx—^f^yY — ^cz^ + {aa+j3b—2yJa(ia)xy , 

(31) 

hence V < Q provided that /36 < 2\/af5a — aa (since 
this would imply —aax^ — + {aa + fih)xy < 0). Let 

us pick /? = aa. Then, the condition reduces to & < 1. 
As a > 0 was arbitrary, we set a = ^. This means that 
whenever a>0,0<6<l, and c > 0 , there exists a 
Lyapunov function 

y{x,y,z) = ]^x^ +'^y^ + , (32) 

since V(0, 0,0) = 0, \V\ —>■ oo as \{x,y,z)\ —>■ oo (ra¬ 
dially unbounded), and V < 0 for {x,y,z) ^ (0,0,0). 
Interestingly, the condition < & < 1 is exactly the sta¬ 
bility condition known in the literature ESI- Therefore, 
parameters implying the existence of this contraction 
map correspond to known stable parameters. 

Now, if we were to seek such a map for only one of 
the state variables, then using what we have obtained 
in Section 3, we find that there exists a contraction map 

Vr/ • -x 1 2 a / 

V{x,x,x) = -x^ + -lx + - j 

^ 2 (33) 

a6 A ai -b (1 -b a)x + ax\ 

2 \ ax J 

for the state variable x{t). Then, one may verify V < 
0 away from the equilibrium a; = 0. One can obtain 
similar contraction maps in either of the other two state 
variables. 


5.3 Qi-Chen-Du-Chen-Yuan System 

In order to find a contraction map for the Qi-Chen-Du- 
Chen-Yuan (QCDCY) system, we begin with the bowl 
shaped assumption for a Lyapunov function about the 
equilibrium {x,y,z) = ( 0 , 0 , 0 ), 

V{x,y,z) = ax"^ + Py'^+-fz'^ , (34) 

where a > 0 , ,5 > 0 , and 7 > 0 are constant parame¬ 
ters to be selected. Differentiating with respect to t and 
using the three constituent equations of the QCDCY 
system, we have 

V = {a — P + 'y)xyz + {aa + Pb)xy — aaj? — Py^ — "fcz^ . 

(35) 

Since we need a > 0, /? > 0, and 7 > 0, we should 
consider model parameters satisfying a > 0 and c > 0 . 


To remove the first term, which is hyperbolic in nature, 
we should choose P = a -b 7 . Meanwhile, the remove 
the second term, which is also hyperbolic, we should 
set P = —^a for non-zero b. Since all other parameters 
are positive, we must require 6 < 0. Then, P = 

and placing this into /3 = a -b 7 gives 7 = . As 

we need 7 > 0 , this gives the added restriction a > | 6 |. 
The parameter a is arbitrary, so we take a = ■^. We 
therefore obtain 




Q- 1^1 ,2 

m ’ 


(36) 


and this candidate function is indeed a contraction map 
satisfying 14(0,0,0) = 0, V < 0 for {x,y,z) ^ (0,0,0), 
and |I4| —>■ 00 as \{x,y,z)\ —>■ 00 , provided that the 
parameter restrictions a > | 6 |, 6 < 0 , and c > 0 hold. 
Therefore, under these parameter restrictions, the zero 
equilibrium is asymptotically stable for the QCDCY 
system. 

When we obtain a single equation for a state vari¬ 
able, even one containing integrals, we may similarly 
obtain a contraction map. Since we have obtained an 
equation for the state variable z{t) in the QCDCY sys¬ 
tem in Section 3, we shall choose to construct a con¬ 
traction map for that state variable here. Doing so, we 
find that 


V{z,z) = 


e2at{^ -b CZp 


1^1 Xq -b 2 /g e^“®(a -b z(s))(i(s) -b cz{s))ds 
-b 2 J e‘^‘^^{a +z{s)){z{s) + cz{s))ds 

a - \b\ 2 

-z 


\b\ 


(37) 


satisfies Y < 0 for all z ^ 0, given that a > | 6 |, 6 < 0 , 
and c > 0. Hence, V{z,z) is a contraction map for the 
state variable z{t) when a > | 6 |, 5 < 0 , and c > 0 . 
With this, we have determined the stability of the zero 
equilibrium for the QCDCY system. 


6 Computational considerations for chaotic 
trajectories 

There are a variety of methods available for trying to 
find chaotic trajectories in nonlinear dynamical sys¬ 
tems, and the approach highlighted in this paper does 
not add to collection of tools, explicitly. However, the 
reduction of dimension approach outlined in Section 2 
can be used to make hnding chaos in dynamical systems 
more efficient. To demonstrate this, we shall consider 









Reduction of dimension for nonlinear dynamical systems 


11 


two rather distinct approaches, namely, the undeter¬ 
mined coefficients method for obtaining chaotic trajec¬ 
tories and the competitive modes analysis for identifi¬ 
cation of chaotic parameter regimes. For each of these 
approaches, we show that an application of reduction 
of dimension results in a simplification of each test for 
chaos. 


6.1 Calculation of trajectories via undetermined 
coefficients 

When attempting to analytically calculate chaotic tra¬ 
jectories, even in an approximate sense, one often re¬ 
duces the dimension of the governing equations. The 
reason for this lies in the fact that it is easier to con¬ 
sider an expansion for one state variable, rather than 
multiple state variables. To best illustrate this point, 
let us return to the Rdssler equation ([3]). 

One popular method for approximating trajectories 
of chaotic systems analytically is the undetermined co¬ 
efficient method [3l|4l[3|6]. Since Taylor series expan¬ 
sions for nonlinear systems often have a finite region 
of convergence centered at the origin, yet the chaotic 
dynamics remain bounded in space, one often consid¬ 
ers non-polynomial base functions. One popular choice 
would be a function of the form 


a)Bn and Cn = —{p?r? -|- aa\n\ -\- l)Bn for all n G 
Z. If we were to truncate the expansion for y(t), in 
the manner described above, we would need to solve 
2 J -|- 1 nonlinear algebraic equations for and 

a, while the coefficients for x{t) and z(t) are immedi¬ 
ately found once we know these parameters. This means 
that by first reducing the dimension of the ODE sys¬ 
tem, we would be able to reduce the computational 
complexity of the problem by a factor of three. For 
higher-dimensional system, the reduction in computa¬ 
tional complexity will scale as the dimension of the sys¬ 
tem itself. In other words, a solution in term of the un¬ 
determined coefficient method will not depend on the 
size of the dynamical system provided that the dynam¬ 
ical system can be reduced in dimension to a single 
equation governing one state variable. 


6.2 Competitive modes analysis: A check for chaos 

The method of competitive modes involves recasting 
a dynamical system as a coupled system of oscillators 
[HIIIIIIIIIMITIIIH]. Consider the general nonlinear au¬ 
tonomous system of dimension n given by 

Xi = fi{xi,X2,.:,Xn) ■ (39) 


Sit; {A,} 


j=oo 
j=-oo ’ 


a) 


fort> 0 , 

forf< 0 . 

(38) 


In this expression, the Aj G R and the parameter a > 0 
are undetermined parameters which one typically will 
obtain in an iterative manner. Assuming such an expan¬ 
sion in time, it males sense to consider a solution the the 
Rossler system Q of the form x{t) = S{t; , a), 

y{t) = S{t- , a), and z{t) = S{t; , a). 

Placing these equations into one would obtain an 
infinite system of nonlinear algebraic equations for all 
of the coefficients and the temporal scaling a > 0. In 
practice, one would truncate these expansions, taking 
the sum over — J < j < J for some J > 0. As the 
solutions may converge slowly - if they converge at all 
(owing to the nonlinearity), one would need to solve 
6 J -b 1 nonlinear algebraic equations. 

Assume, instead, that we wish to solve o by the 
approach described above. We would then insert the 
expansion for y{t) into ([4]). Assuming that we can solve 
the resulting nonlinear algebraic equations for the con¬ 
stants and a, we can then recover x(t) and 

z{t) by recalling x = y — ay and z = —y + ay — y. From 
these expressions, it is simple to show A„ = — (a|n| -I- 


Differentiation of (IMl) once gives a coupled system of 
second order equations. 


Xi — 


n 


E/. 


dxj 


T (xi, X2 , ■. ■, ^i—1 , , ..., Xn} . 

(40) 


When a gi is positive, its respective zth equation be¬ 
haves like an oscillator. The following conjecture is posed 
in [H] : 

Competitive Modes Requirements: The conditions 
for dynamical systems to be chaotic are given by: 

(A) there exist at least two modes, labeled gi in the 
system; 

(B) at least two g's are competitive or nearly competi¬ 
tive, that is, for some i and j, gt ^ gj > 0 at some t; 

(C) at least one of the g’s is a function of evolution 
variables such as t; and 

(D) at least one of the h’s is a function of system vari¬ 
ables. 


The requirements (A)-(D) essentially tell us that a 
condition for chaos is that two or more equations in (1401) 
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behave as oscillators {gi > 0), and that two of these os¬ 
cillators lock frequencies at one or more times. In prac¬ 
tice, we find that the frequencies agree at a countably 
infinite collection of time values [1^118] . The frequen¬ 
cies should be functions of time (i.e., we have nonlinear 
frequencies), and there should be at least one forcing 
function which depends on a state variable. 

In order to consider all possible chaotic dynamics, 
one would have to compare each pair gi = gj, i j, 
i,j = l,2,3,...,n. Accounting for symmetry, this gives 
2”“^ — 1 matchings to consider. For high-dimensional 
dynamical systems, this number becomes rather large. 
As an example, in the case of a ten-dimensional system, 
there will be 511 possible matchings to be checked in 
order to ensure one has determined the possible chaotic 
regimes. For such situations, the approach is not par¬ 
ticularly efficient, and a competitive modes analysis is 
often considered for systems of dimension three or four 
in the literature. 

Let us consider dynamical systems dMl) which can 
be put into the form of a single equation for one state 
variable. As such an equation will encode the dynam¬ 
ics of the complete system, it is sufficient to consider 
a competitive modes analysis for the resulting equa¬ 
tion. Suppose that the resulting equation has maximal 
derivative of order p > 0 (where p need not be equal to 
n, as we have seen in earlier sections). Then, associating 
y{t) to this single state variable, we have 


dtP ~ y' dt''"' dtP-^ ) 


(41) 


Since the competitive modes analysis relies on us ob¬ 
taining a system of oscillator equations, let take y = yi 
and write the equation (HD) as the system 


yi = y 2 ■ • ■, ijp-i = yp,yp = A(yi, ...,yp). (42) 

Differentiation of (HU) results in the system of second 
order equations given by 


i/1 = 2/3, 


i/p—2 — 2/p) 

i/p-1 =yp = F{yi,...,yp), 


.. ^ OF dF ^. 

yp = 2^ -^yi = 2^ —2/*+i + —^(//v • ■ ■, i/p) • 






dyF 


dyp 


(43) 


The right hand side of the first p—2 second order equa¬ 
tions do not depend on the state variable for each re¬ 
spective equation, so 51 = • • • = gp -2 = 0. Hence, these 


equations are never oscillators. Meanwhile, we can de¬ 
compose the right hand sides of the latter two equa¬ 
tions, so that 

Fiyi, ■■■,yp) = -2/p-iSp-i + hp-i (44) 

and 

dF 

F ^2/*-si + ■■■^yp) = -yp9p + K ■ (45) 

Note that there are now exactly two mode frequen¬ 
cies, gp-i and gp, and there is always an hi depend¬ 
ing on state variables. In order to determine if the sys¬ 
tem satisfies the competitiveness conditions (A)- 
(D) (and therefore if the original system satisfies these 
competitiveness conditions), it is sufficient to check if 
gp-i = gp> 0 for some collection of time values. This is 
only one condition to check, rather than 2"“^ — 1 con¬ 
ditions to check from the original system. Therefore, 
while the conversion of the system (l39ll to the equiva¬ 
lent system (l42ll may seem somewhat roundabout, do¬ 
ing so greatly simplifies the search for possible chaotic 
dynamics under the competitive modes framework. 


7 Discussion 

The construction of Lyapunov functions for nonlinear 
dynamical systems is often either simple, or quite chal¬ 
lenging, with little room in between. Aside from choos¬ 
ing some standard forms (such as the common bowl 
shape centered about an equilibrium value), there is 
more an art to the selection of such function. However, 
as we have demonstrated for the Rossler system, it is 
possible to use differential elimination to obtain a con¬ 
traction map in a single state variable, which can then 
be used to obtain a stability result for all state vari¬ 
ables. A similar approach was also employed in [1] to 
study Michaelis-Menten enzymatic reactions, and after 
using a reduced equation, the proof of global asymp¬ 
totic stability for a positive steady state was rather 
simple. Differential elimination, and reduction of di¬ 
mension more generally, can therefore be useful in help¬ 
ing one determine asymptotic properties of solutions to 
nonlinear dynamical systems. On the other hand, it is 
now known that there are autonomous chaotic systems 
which lack any equilibrium [78] . The approach outlined 
here could still be used to construct maps which demon¬ 
strate the boundedness of trajectories in time for such 
systems, even if such trajectories do not approach any 
fixed points. 

While useful for obtaining contraction maps, reduc¬ 
tion of dimension is also a promising tool for finding 
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and studying chaotic trajectories in nonlinear dynam¬ 
ical systems. We were able to show that reduction of 
dimension can be used to simplify calculations involved 
in using undetermined coefficient methods [3l|4l[5l[6] by 
a factor of 1/n, where n is the dimension of the dy¬ 
namical system, meaning that the number of equations 
needing to be solved will not increase with the size of 
the system - i.e. the computational complexity of the 
approach will not scale with the size of the system but 
rather will remain fixed. This means that one may ap¬ 
proximate chaotic trajectories through such approaches 
in systems of rather large dimension, without the com¬ 
putational problem becoming unwieldy, provided that 
the system can be reduced to a single ODE governing 
only one state variable. 

While chaos is often studied numerically, recently 
analytical approaches have been employed to construct 
trajectories approximating chaotic orbits. When em¬ 
ploying these methods, it can be beneficial to consider 
a single equation rather than a system of equations, 
even if the single equation is more complicated. This is 
true for series and perturbation approaches, as the re¬ 
duction of order requires one to track fewer functions, 
which is particularly useful when dealing with messy 
equations. Furthermore, analytical approaches permit¬ 
ting the control of error, such as the optimal homotopy 
analysis method, rely on assigning an error control pa¬ 
rameter to each state variable. Reduction of dimension 
can allow one to minimize error via a single control pa¬ 
rameter [5], rather than over multiple parameters ca, 
which is computationally less demanding. 

The reduction of dimension can also be useful for di¬ 
agnostic tests for chaos. When applying the competitive 
modes analysis, a type of diagnostic criteria for finding 
chaotic in nonlinear dynamical systems, one performs 
binary comparisons between the mode frequencies of 
an oscillator corresponding to each equation. However, 
if one first applies reduction of dimension and reduces 
the system of a single ODE for one state variable, we 
prove that only one comparison would needed. This is 
particularly beneficial when studying large dimensional 
systems, since the number of naive comparisons needed 
scales like 2"“^ in dimension n. 

In the future, it would be interesting to consider 
an algorithm that considered elimination not only ele¬ 
ments 9" (9 = Jj) of a differential ring but also 

integral operators (where 9“" satisfies = 9° 

and hence is the inversion of the operator 5"). Indeed, 
in cases where the differential elimination algorithm 
failed to give a reduction of the system to a single ODE, 
we found by manual substitutions that one can arrive 
at an integro-differential equation. While more compli¬ 
cated, such integro-differential equations can still cast 


light on the behavior of solutions, and can prove useful 
in obtaining Lyapunov functions. More generally than 
for dynamical systems, these inverse operators play 
a role in the study of operators and integrable hierar¬ 
chies arising in nonlinear evolution PDFs [SOllSIllMlIS^ 
[MII85U86) . Therefore, the extension of the algorithm to 
the ring of formal Laurent series in d would be a fruit¬ 
ful area for future work, not only for dynamical systems 
but also for integrable partial differential equations. 
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Appendix A: List of Chaotic Systems 

When testing our approach, we considered a variety of specific 
chaotic systems. Results for these are listed in Table 1. As the 
form and scaling of such equations can vary in the literature, 
we list the specific form of these equations used in our work. 

Let x,y,z,w G C"(R) where n is the dimension of the 
relevant dynamical system, and let a,b,c,d G K be parame¬ 
ters. 

Lorenz system ElESI: 


x = a{y -x) , 
y = x{b-z)-y, 
z = xy — cz . 

Modified Chua’s circuit f? l?ni38l : 

X = a — (2x^ — 

y = X - y + z , 
z = -by. 

Chen - Lee system IMlIiOl : 

X = ax — yz , 
y = by+ XZ, 

1 

z = cz —xy . 

3 

Rabinovich - Fabrikant equations 
X = y (^z — 1 -\- x^') -\- ax , 
y = X (Zz-\-1 - x^) -\- ay , 
z = —2z{b + xy) . 

Rossler system [HES]: 
x = -y - z, 
y = x + ay, 
z = b + z{x — c) . 


(46) 


(47) 


(48) 


(49) 


(50) 
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Chen system SSISZ]: 


Appendix B: List of Hyperchaotic Systems 


X = a{y - x) , 

y = (b — a)x — xz + by , (51) 

z = xy — cz . 

Lii system I48II49I : 

X = a{y - x) , 

y = by-xz, (52) 

z = xy — cz . 

T system I5f)linil52l : 

X = a{y - x) , 

y = {b — a)x — axz , (53) 

z = xy — CZ , 

4D Qi-Du-Chen-Chen-Yuan system m- 

X = a{y — x) + yzw , 
y = b(x + v) — xzw , 

^ ^ (54) 

z = —cz + xyw , 

w = —dw + xyz . 

Qi-Chen-Du-Chen-Yuan system 1541 : 

X = a{y - x) + yz , 

y = bx — y — xz , (55) 

z = xy — cz . 


Generalized Lorenz canonical form l55l[56l[57l : 

X = ax — [x — y)z , 

y =-by - {x - y)z , (56) 

i = -cz + (a; + y){x + dy) . 


Two-parameter model for the blue-sky catastrophe 1581 

mm- 

1 = (2 -h a - 10 -h j/^)) Xy'^ + 2y z'^ , 

y = - (1 + 2/) {y^ +‘^y + -4:X + ay, (57) 

2 = (1 + y)z^ -h - 6. 

4D Lorenz - Stenflo system l6T1l62ll63l : 


X = a{y — x) + bw , 
■y = cx — xz — y, 
z = xy — dz , 
w = —X — aw . 


(58) 


Genesio-Tesi system 15111551 : 
x = y, 

y = z, (59) 

z = ax + by + cz + x^ . 

Arneodo-Coullet-Tresser system [55]: 

x = y, 

y = z, 

z = ax — by — cz — x^ . 


When testing our approach, we considered a variety of specific 
hyperchaotic systems. Results for these are listed in Table 2. 
As the form and scaling of such equations can vary in the 
literature, we list the specific form of these equations used in 
our work. 

Let X, y, z,w E C^(R) and let a, b, c, d,e, f £ R be param¬ 
eters. 

4D Rossler flow 1671 : 


x = -y- z, 
y = X 0.2hy -R w , 
z = S + xz , 
w = —0.52 -|- 0.05ui. 

Hyperchaotic Chen system [BSlil]: 

x = a{y-x), 
y = —bx — xz + cy — w , 
z = xy — dz , 
w = X . 


(61) 


(62) 


Hyperchaotic Lii system ED): 

X = a{y — x) -|- ui, 
y = by- xz, 
z = xy — cz , 
w = xz + dw . 

Modified hyperchaotic Lii system US- 

X = a(y - X + yz) , 
y = by — xz + w , 
z = xy — cz , 
w = —dx. 

Hyperchaotic Wang-Liu system [72]: 

x = a{y-x), 
y = bx — cxz -|- w , 
z = —dz -|- ex^ , 
w = —fx. 

Hyperchaotic Jia system Em]: 

X = a(y — x) -|- ui, 
y = bx — xz — y , 
z = xy — cz , 
w = dw — xz . 


(63) 


(64) 


(65) 


( 66 ) 


Hyperchaotic Qi - van Wyk - van Wyk - Chen system 1751 

E5]: 


X = a{y - x) + yz , 
y = b{x + y) - xz , 
z = —cz — dw -|- xy , 
w = cz — fw -|- xy . 


(67) 


(60) 












